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Abstract. We present a method to compute the degrees of the Segre classes 
of a subscheme of complex projective space. The method is based on generic 
residuation and intersection theory. It has been implemented using the soft- 
ware system Macaulay2. 



1. Introduction 

Segre classes are generalizations of characteristic classes of vector bundles and 
they occur frequently in intersection theory. Many problems in enumerative geome- 
try may be solved by computing the Segre classes of an algebraic scheme. Given an 
n-dimensional subscheme Z of complex projective space there are n + 1 Segre classes 
of Z. The i th Segre class is a rational equivalence class of codimension i cycles on Z. 
Thus a Segre class may be represented as a weighted sum of irreducible sub-varieties 
Vi, . . . , V m of Z. The degree of a Segre class is the corresponding weighted sum of 
the degrees of the projective varieties V\, . . . , V m . 

In this paper we present a method to compute the degrees of the Segre classes 
of Z , given an ideal defining Z . The procedure is based on the intersection theory 
of Fulton and MacPherson. More specifically, we prove a Bezout like theorem 
that involves the Segre classes of Z and a residual scheme to Z. The degree of 
the residual may be computed providing intersection-theoretic information on the 
Segre classes. This enables us to compute the degrees of these classes. 

If Z is smooth, the degrees of the Segre classes of Z carry the same information 
as the degrees of the Chern classes of the tangent bundle of Z. For instance, when 
Z is smooth, the degree of the top Chern class, which is equal to the topological 
Euler characteristic of Z, can be computed from the degrees of the Segre classes of 
Z. Hence one can compute the topological Euler characteristic using our procedure 
(provided that the input ideal defines a smooth scheme) . The relationship between 
the degrees of the Segre classes of Z and the Chern classes of the tangent bundle of 
Z may be extended to the non-smooth case via the so-called Chern-Fulton classes 
of Z. In case Z is smooth, the Chern-Fulton classes coincide with the Chern classes 
of the tangent bundle. 

We would like to mention two interesting features of our method. One feature 
is that the algorithm is completely elementary in that it requires no knowledge of 
intersection theory to understand the steps in the procedure and only very basic 
background in algebraic geometry. The method therefore provides a way of under- 
standing the computation of Segre classes from an elementary point of view. In 
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addition, the algorithm is easy to implement. Another feature is the fact that our 
method is implcmentable in a numerical setting via numerical homotopy methods, 
see [TT] for an overview of this area. This allows the method to be applied in settings 
that can be time consuming (and even out of reach) of current symbolic methods. 
One such setting is when the generating set, for an ideal determining Z, has com- 
plicated coefficients. An additional setting where current numerical methods can 
sometimes obtain useful information about the degrees of Segre classes, beyond the 
reach of current symbolic methods, is when Z is a reduced scheme of high codimen- 
sion. These gains come through numerical approximation and parallclization (but 
at the expense of exactness). 

The procedure presented in this paper has been implemented in the sym- 
bolic setting using the software system Macaulay2 7 and is available at 
http:/ /www. math. su.se/~jost/segreimplementation. htm. A numerical version, uti- 
lizing the software package Bertini [3], will be made available pending the release 
of Bertini Version 1.3. Initial experiments, involving subvarieties of relatively large 
dimension and codimension, show a great deal of promise for the algorithm in the 
numerical setting. 

In the paper [I] Alum formulates an algorithm that also computes the degrees of 
the Segre classes of a subscheme Z of projective space. In addition he shows how 
to reduce the computation of the so-called Chern-Schwartz-MacPherson classes of 
a subscheme of projective space to the computation of the degrees of certain Segre 
classes. In the present paper we present an alternative method to Alum's. Though 
the two are closely related, they have a rather different computational behavior 
and seem to complement each other well (see Section [6] for some examples) . Apart 
from the difference in computing speed in various cases one may ask what is the 
need for another method with the same output as an existing method. One answer 
is that our approach is different and therefore sheds new light on the problem of 
computing Segre classes. But more importantly we would answer by repeating the 
two features mentioned above, namely that our method is elementary and that it 
is readily amenable to numerical computation. 

The paper is organized as follows. In Section [2] we give the basic definitions and 
state a theorem from intersection theory. In Section[3]we derive a recursive formula 
for Segre classes which is the basis of our method. The procedure to compute Segre 
classes is presented in Section |U Some examples are given in Section [5] and in 
Section [5] we give a list of run times on examples comparing our method to other 
algorithms. 

The results of this paper are generalizations and variants of the results in 2,4 
to the setting of subschemes of projective space. 
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2. Background in intersection theory 

We start by going through some concepts and results from intersection theory. 
For this paper, the main reference on matters of intersection theory is Fulton's book 
0. 
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2.1. Notation. Let Y be an algebraic scheme over C of dimension n. By a sub- 
scheme of Y we will mean a closed subscheme. We will denote by C P (Y) the free 
Abelian group on irreducible p-dimensional subvarictics of Y. The p th Chow group 
of Y is the quotient of C p (Y) by the cycles rationally equivalent to 0, and it is 
denoted A p (Y). The Chow group of Y is the group A*(Y) = 0™ =o A p (Y). Given 
an element a 6 A*(Y), {a} p € A p (Y) will denote the p th homogeneous compo- 
nent of a (if n < p, then {a} p — 0). A subscheme X C y induces a cycle class 
[X] 6 A*(y) represented by X)*=i TO i^, where Xx, • • . ,X 4 are the irreducible com- 
ponents of X and m\, . . . ,m t their geometric multiplicities in X. In particular, 
[0] =0. If a e A* (X), we will at times consider a to be an element of A*(Y), 
omitting in the notation the push-forward under the inclusion map. 

For a rank p vector bundle E on Y we have the Chern class operations Ci(E) : 
A p (Y) -> A p -i(Y) for i < p, < i < p and < p < n, see [5] Chapter 3. 
The value of Ci(E) on a € -A P (Y) is denoted ci{E) n a. The corresponding map 
A*(Y) A*(Y) is also denoted c^E 1 ), where Ci{E)C\a — if a e and p < i. 

The total Chern class operation c(£') : A*(y) — > A*(y) : a i-» c(E)C\a is defined by 
c(E) = J2i Ci{E). The operation of a product of Chern classes on the Chow group is 
defined as the composition of the individual Chern class operations. The map Cq(E) 
is the identity homomorphism. If Y is smooth, the Chern classes are well defined 
elements of A*(Y) and the operations Ci(E) n a correspond to the intersection 
product. Let X C Y be a closed subscheme with inclusion i : X — > Y and let 
a G A*(X). By the notational convention mentioned above we will sometimes 
write c(E) Ha to mean c(E) n i*{a) £ A*(Y). By the projection formula, see [5] 
Theorem 3.2 (c), c(E) n i»(a) = i*(c(i*E) n a). 

For a Cartier divisor D on y, the corresponding line bundle on y is denoted by 
O y {D). 

In this paper, varieties are by definition irreducible and reduced. Finally, we use 
the convention dim(0) = — 1. 

2.2. Regular embeddings. Let Y be an algebraic scheme over C. A closed em- 
bedding X —> Y of a subscheme X of Y is called a regular embedding of codimension 
d if the following holds. Every point of X has an affine open neighborhood U in Y 
such that the ideal defining X n U is generated by a regular sequence of length d 
in the coordinate ring of U. 

Lemma 2.1. Let Y be a complex variety and let C be a line bundle on Y such 
that the corresponding complete linear system is base point free. For global sections 
o~i, . . . , a p £ r(y, C) of C, let X = X(o~i, . . . , a^) denote the scheme of common 
zeros of ctx, ■ • ■ , c P - Then, for general o~\, . . . , S Y(Y,C), the natural embedding 
X — > y is regular. 

Proof. Let <f> : Y — > W be the map given by the complete linear system correspond- 
ing to C and put e = dim(im(</>)). If jj, > e then A(t7x, . . . , a^) = for general 
sections a\,. . . ,a p and hence we may assume that p < e. That X(o~i) — > Y is 
a regular embedding for a general section g\ is clear since Oy(U) is an integral 
domain for any open set U C Y. Hence, for any affine open U C Y, Oy(U) has 
no zero divisors and any ui / will give a regular sequence in Oy(U). If p = 1 
we are done. If dim(Y) < 1 then e < 1 and therefore = 1 in this case. Assume 
that 1 < p, in particular 1 < dim(y). Then e > 2 and it follows from Bertini type 
theorems that X(a%) is a variety for a general section ox. In fact, A(crx) is reduced 
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by [6] Corollary 3.4.9 and X(a\) is irreducible by [6] 3.4.10. Replacing Y by X(ai) 
and restricting C to X{a\) we have reduced to the case of a lower dimensional 
ambient variety since the composition of two regular embeddings is regular (see [3] 
Appendix B.7.4). □ 

2.3. Segre classes and intersection products. Let Y be a complex variety and 
let X be an n-dimensional subscheme of Y. Suppose that 1^7 and let Y be the 
blow-up of Y along X. Let tt : Y —> Y be the projection, let X = tt^ 1 (X) be the 
exceptional divisor, and let r\ = 7r|jf . The total Segre class s(X, Y) of X in Y is an 
element of A*(X) which may be characterized as follows (see [5] Corollary 4.2.2): 

s (i,y) = ^(-irV(i p ). 

P>1 

Here X p is the self intersection of Cartier divisors defined in [5] Definition 2.4.2. 

Remark 2.2. The Cartier divisor X on Y corresponds to a line bundle on Y whose 
restriction to X is the normal bundle of X in Y. The dual bundle to the normal 
bundle is denoted 0(1). The exceptional divisor X is naturally identified with the 
so-called projective normal cone P{CxY) of X in Y. The total Segre class is given 

by 

S (X i y)=^ry t ( Cl (0(l)rn[P(CxF)]). 

i>0 

This definition of Segre classes generalizes beyond normal cones of subschemes to 
arbitrary cones, see [5] Chapter 4. 

We will now recall the Fulton-MacPherson approach to intersection products, 
see [5] Chapter 6. Let Y and V be complex varieties and put k = dim(V). Let X 
be a closed subscheme of Y such that there is a regular embedding i : X — > Y of 
codimcnsion d < k. Let / : V — > Y be a morphism and put W = f~ 1 (X). Then 
we get the following fibre product diagram 

w— ] -^v 

9 / 

X r^Y 

i 

where j : W — > V is the inclusion and g : W — » X the restriction of / to W . The 
normal cone CwV of 14 7 in V can be constructed as follows. Suppose first that V 
is affine with coordinate ring A and that W is defined by an ideal J C A generated 
by /i, . . . , fd £ A. Let B = A/ J. Then CwV is the spectrum of the B-algebra 
©s>o J s I J s+1 ■ Thus CwV may be embedded as a closed subscheme of W x C d 
defined by the kernel of the surjective homomorphism 

B[xi s ...,x d ]->0JVJ' +1 

s>0 

which maps a;, to the image of /, in J/ J 2 . In the general case, CwV may be 
constructed by covering V with open affine subsets and gluing the normal cones 
of the affine patches together. In case the embedding W — > V is regular, CwV 
is a vector bundle, namely the normal bundle NwV. In this case the total Segre 
class of W in V is the inverse of the total Chern class of NwV in the sense that 
s(W, V) = ciNwVy 1 n [W], where c(NwV)- 1 is a formal inverse of c(N w V). The 
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normal cone CwV has pure dimension k, see [5] Appendix B.6.6. Let C = CwV, 
let N = g*(N x Y) (where N X Y is the normal bundle of X in Y), and let p : N -> W 
be the projection. It is shown in [5] Chapter 6, that C embeds in N and therefore 
it determines a class [C] € Ak(N). Now, by [5] Theorem 3.3 (a), the flat-pullback 
p* : Ak-d{W) — > Ak(N) is an isomorphism. This map is given by p*(Z) = [p^ 1 (Z)] 
for an irreducible subvariety Z C W. The intersection product of V by X on Y is 
a class in Ak-d(W) denoted X ■ V and defined by 

x.v=(j>*r 1 (m 

An important connection to Segre classes is given by [5] Proposition 6.1 (a): 

X-V = {c(N)na(W,V)} k - i . 

We will now state the Residual Intersection Formula from [5] which is the main 
result underpinning our method to compute Segre classes. Let Y, X, W , V, N, k 
and d be as in the above definition of the intersection product X ■ V. Let Z C W 
be a closed subscheme and suppose that Z ^ V. Let ir : V — > V be the blow-up 
of V along Z and put W = 7r -1 (W) and Z = 7r _1 (Z). Let be the residual 
scheme toZinW with respect to F, see [5] Definition 9.2.1. This is a scheme such 
that, if ,y{Z), J^(W) and ^{R) are the ideal sheaves in 0^ defining the respective 
schemes, then 

J r (W) = J r (Z)-J r (R). 
Let r\ : W — > W be the restriction of 7r to W and let O(-Z) denote the pullback of 
Oy(— Z) under the inclusion W — > V . 

The following proposition is Corollary 9.2.3 of [5]. 

Proposition 2.3. With notation as above, 

X-V = {c(N)ns{Z, V)} k - d + R, 
where K = r)*({c(r]*N ® O(-Z)) n s(5, F)} fe _ d ). 

3. Computing Segre classes of projective schemes 

Let Z be a proper n-dimensional subscheme of complex projective space P . This 
paper is about a method for computing the push-forward of s(Z, P fc ) to P fc given 
an ideal defining Z. In this section we explain how to derive information about the 
push-forward given sufficiently general elements from the ideal. 

Let So, . . . , s n be the homogeneous components of s(Z, P fc ) with Sj of codimension 
i, that is s(Z, P fe ) = Y^i=o Si w ^ tn Si e ^n-i(Z)- The degree of a 0-cycle a = 
J2i ^Pi on P fe , m, € Z and pj g P fe , is simply deg(a) = ^ m*. If 7 : Z -> is 
the inclusion map, we define the degree of Sj by 

deg( Sl ) = deg( 7 *(^)-ff M ), 

where H € Afc_i(P fe ) is the hyperplane class and the product is the intersection 
product on P fc . The numbers {deg(si)}i is the output of our procedure, they carry 
the same information as the push-forward y*(s(Z, P fc )). The degree of any a £ 
A p (P k ) is defined similarly by deg(a • H p ). 

Let I C C[xq, ■ ■ ■ , Xfc] be a homogeneous ideal. For a positive integer m, we use 
I(m) to denote the m th graded piece of /. Given a homogeneous ideal J, the ideal 
quotient J : I is given by 

J:I = {f£C[x ,...,x k ]:fICJ}, 
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and the saturation of J with respect to 7 is 

J : I°° = (J J : P. 

p>i 

Note that (J : P) : I = J : P +1 for p > 1, that the ascending sequence of ideal 
quotients J : I C J : I 2 C J : I 3 C . . . stabilizes and that J : P = J : I°° for large 
enough p. 

Remark 3.1. Let J and J be homogeneous ideals of C[xq, . . . , Xk] and let V(J) 
and V(J) denote the corresponding zero-loci in P fe . If I = C[xo, ■ ■ ■ , Xk], then 
J : 7°° = J. Suppose 7 ^ C[xo, . • . , x&]. The ideal J : J°° is homogeneous and 
the scheme R defined by J : 7°° is supported on the Zariski-closure of V(J) \ V(I). 
In fact, if J = Pi; Qi is a primary decomposition (so each Qi is homogeneous and 
primary) then 

(1) J:7°°= f) Qi. 

{i:V(Q z )£V(I)} 

To see this, note that J : 7°° = Hi(Q« : /0 °) and that £ ^(-f) precisely when 

3 7. If V'Qi" 3 7, then Qj D 7 P for some p and Qi : 7°° = C[xo, • ■ • , Xk}- On 
the other hand, if y/Qi 7, then Qi : I = Qi since Qi is primary. It follows that 
Qi : 7°° = Qi in this case. 

The following theorem is a Bezout like equality which gives rise to a recursive 
formula for the degrees of the Segre classes of Z in P k . Using the statement of 
the theorem, we may express the degree of a Segre class s p in terms of deg(sj) for 
i < p and the degree of a certain residual scheme R. Computing the degree of the 
residual R is the main computational step in the method. 

Theorem 3.2. Let Z C P k be a subscheme of dimension n defined by a non-zero 
homogeneous ideal I C C[xo, . . . , Xk] ■ Let so, . . . , s n denote the Segre classes of Z 
in ¥ k . Let go,...,g r be a set of non-zero homogeneous generators of 7 and put 
m = maxi{deg(gi)}. For k — n < d < k and general elements fx, ■ ■ ■ , fd G 7(m), 
the following holds. If J is the ideal generated by {/ l5 . . . , /^} and R C P fe is the 
subscheme defined by J : I°° , then 

m d = deg(i?)+^ _.)m*'- i deg(« i ), 
i=o \P l ' 

where p = d — (k — n). 

Proof. The proof is divided in steps through 4. 

Step 0: setup. Let 7' be the ideal generated by 7(m) and let m = (xo, ... ,Xk). 
Then 7 : m°° = 7' : m°°, and therefore 7 and 7' define the same scheme Z CP k . 
We may thus assume that go, . . . ,g r all have degree m. Let 

7T : P fc -> P fc 

be the blow-up of P k along Z and put Z = ^(Z). The map P k \ Z -> P r defined 
by go, ... , g r extends to a map 

cj) : P k -+ P r , 

see 4.4. In fact, P k embeds in P k x P r in such a way that (P fc \ Z) is the graph of 
the map P k \ Z — > ¥ r and <j> is the projection. Let W C P fc be the scheme defined 
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by /i, . . . , fd and put W = ir 1 (W). Let R be the residual toZinW with respect 
to P fc . 

Step 1: we will show that R — »• F k is a regular embedding and that R is ei- 
ther empty or of pure dimension k — d and that no irreducible component of R is 
contained in Z. By [5] 4.4., 

= 7T*((IV(m))®0p fc (-Z). 

In concrete terms, P fe is defined by a bi-homogeneous ideal 

If c Cfio, . . .,xk][yo, -..,y r ] 

such that K contains the elements gijjj — gji/i for < i < j < r. Observe that Z is 
given by the vanishing of go , . . . , g r . Consider the affine open set 

U = U al3 = {(x , ...,x k ,y ,...,y r )eF k :x a ^ O,y ^ 0} 
and let wo = ■ ■ ■ iW r = ^ with wp = 1. Then wq, . . . , wp, . . . , w r are coordi- 
nates on C = {j/^j 7^ 0} C P r . Note that for all i, fi = J2]=o f° r a general 
vector (A°, . . . , A[) € C +1 . With an abuse of notation, we use gt, fi and Wi to de- 
note the corresponding elements of the coordinate ring of U . Then, gj = Wjgp for 
all j. Hence ZC\U is defined by gp. Also, fi — (J2j=o H w j)9p- It follows that RnU 
is defined by the ideal (X^i=o A-ju;., , . . . , X)j=o ^d w i)- We conclude that i? = ijf> _1 (-L) 
for a general linear subspace L C P r of codimension d (if r < d, then i? = 0). Hence 
i? is either empty or of pure dimension k — d and RP\ Z is either empty or of pure 
dimension k — d — 1. It follows that no irreducible component of R is contained 
in Z. Since P fc is a variety (see [5] Proposition II. 7. 16 or [5] Appendix B.6.4), it 
follows by Lemma HOI that the embedding of R in F k is regular. 

Step 2: applying Proposition 12.31 Let X\,..., Xd be defined by X v = {/„ = 0}. 
Then W = fllLi x f ■ Let x = x i x ' ' ' x x d and Y = F k x • • • x P fe (d factors). Let 
j : W — » P fc be the inclusion and let / : P fe -> Y and g : — >• X be the diagonal 
morphisms. The morphism i : X — > Y induced by the inclusions X%, . . . , Xd C P fc 
is a regular embedding of codimension d. Put N = g*{NxY). Letting V = P fc , we 
apply Proposition 12 . 31 to the diagram 

rLi*» — 

X 1 x • • • x X d — P fc x • • • x P fc 

and conclude that 

(2) X ■ F k = {c(N) n s(Z, F k )} k - d + M. 

Here s(Z, F k ) is regarded as a rational equivalence class on W. 
Step 3: we will show that 



p / . 

m<i = ^ I ■ ) mV ~ l de S(s t ) + deg(j, 



i=0 



p — I 



wherep = d- (k-n). We shall first see that N = j*(E) where E = @ v=1 O f k(m). 
Let p v : Y —> F k be the v projection and let X„ be the divisor p~ 1 (X u ) on Y. 
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Then, by [5] B.7.4, 

d 

n x y = @ev(x l/ )u. 

!/=l 

Since /*((9y(X I/ )) = Opfc(m) for all we have that 

j*(0 ¥k (m)) = (/oj)*(£V(X„)) = (io 9 )*(O r (X I/ )) = 9 *(Or(X,)U) 
for all ^. Hence 

,/ 

N = g*(N x Y) = ®g*(0 Y (X u )\ x ) = j*E. 

Note that c(E) = (1 + mH) d € A*(P fe ), where e A k _ 1 (¥ k ) is the hyperplane 
class. 

We now push both sides of ^ forward to P fe by j and then take degrees. By 
Bezout's theorem, deg(j*(X-P )) = m (see [5] Example 6.2.6). By the projection 
formula, j*(c(N) n s{ZJ k )) = c{E) n (s(Z, P fe )). We get that 

i*({c(iV) n s (z,P fc )} fe _ d ) = {j*(c(iV) n s(Z,F k ))} k _ d = 

n 

{(1 + m.H) d ■ j,{s(Z, P k ))} k - d = {(1 + mH) d ■ £ j,(«0}fc-«I- 

i=0 

The degree of the latter expression is 

^m p ^deg( Si ), 

where p = d — (k — n). 

Step 4: it remains to see that deg(j*(IR)) = deg(i?) where R C P fc is the scheme 
defined by J : I°°. In fact,jwe shall see that = [iZ] in A»(P fc ). Let 77 : W -> W 
be the restriction of 7r to W. Since i? — > ¥ k is a regular embedding we have that 
s(R, P k ) — c(N^P k )~ 1 n [R]. Since R is either empty or has pure dimension k — d, 

{c(rfN ® O(-Z)) n fl (5, P fc )}fe-d = 

{^tv <g> O(-z)) n (c(JV^)- 1 n = [i?]. 

Hence R = i]*([R\). Let . . . , Rt be the irreducible components of R and let 
mi,...,m f denote their geometric multiplicities. Since none of the components 
Ri, . . . , Rt is contained in Z and it : P k — > P fc is an isomorphism outside Z, j*(K) = 
Si=i m A' K {Ri)\- Observe that it induces an isomorphism (R \ Z) = (W \ Z). It 
follows from ((J) that £* =1 W 7 ^)] = [&}■ □ 

Remark 3.3. With notation as in the proof of Theorem 13. 2[ note that the group 
GL(C fc+1 ) x GL(C r+1 ) acts transitively on P k x P r and that P fc is regular outside 
Z. It follows by Kleiman's transversality theorem [10] that for a general linear 
subspacc L C P r of codimension <i, i? = P fe n (P k x L) is regular outside Z and the 
multiplicities mi, . . . , m t of the components of i? are all equal to 1. Moreover, it 
follows that the scheme R defined by J : 7°° is regular outside Z. This could be of 
interest in connection with the computation of the degree of i?, which is the main 
computational ingredient in our method to compute Segre classes. 
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4. The method 

Thcorem l3.2l states that certain conditions hold for a general choice of elements of 
a given ideal. By choosing these elements randomly we turn this into a probabilistic 
algorithm. Applying Theorem 13.21 to solve for the Segre classes recursively, we 
obtain the following procedure to compute the degrees of the Segre classes of a 
subscheme of projective space. The input is an ideal defining the subscheme. 



Procedure 1 A method to compute the degrees of Segre classes 

Input: Non-zero homogeneous generators g$, . . . , g r of an ideal / C C[xo, . . . , Xfe]. 
Output: The degrees of the Segre classes of the subscheme of P fe defined by I. 

1: Let m = maxi{deg(gi)}. 

2: Let Z C P k be the scheme defined by I and compute n = dim(Z). 

3: Pick random elements /i, • • • , /fc G J(m). 

4: for d = k — n to k do 
5: Let J= 

6: Compute deg(i£), where R C P fe is the scheme defined by J : I°° . 
7: Let p = d — (k — n) and compute 

deg( Sp ) = m d ~ deg(R) - W M deg( Sl ). 

8: end for 

9: return deg(s ),..., deg(s„) 



Remark 4.1. We will use the notation of Procedure [T] and the proof of Theorem l3.2l 
In particular tt : P fe — > P fc denotes the blow-up of P fe along Z. Instead of saturating 
with respect to the whole ideal I, as is done in Procedure [1] one could saturate 
with respect to one element of I. Let h € I, h ^ 0, let R' be the scheme defined 
by J : (h)°° and let C P fc be the hypersurface defined by h. The claim is that 
we could replace R by R' in Procedure [T] Tracing back the conditions on R used 
in the proof of Theorem 13.21 we see that we only need to show that no irreducible 
component of the residual R C P fc is contained in 7r -1 (if). This follows exactly as 
in step 1 of the proof of Theorem 13.21 where it is shown that R has no irreducible 
component inside the exceptional divisor Z . 

Remark 4.2. Let Z C P fc be a subscheme of dimension n and let Si € A n ~i(Z) be 
the Segre classes of Z, that is s(Z, P fe ) = X)"=o Si - Define the total Chern-Fulton 
class of Z by 

c'(Z)=c(T pk \ z )ns(Z,¥ k ), 

where T V k is the tangent bundle of P fe . This definition is independent of the embed- 
ding of Z in P fc in the sense that if Z admits two embeddings into smooth varieties 
M and P, then c{T M \z) H s(Z, M) = c{T P \ z ) n s(Z, P), see [5] Example 4.2.6. Let 
c'{Z) = Y^i=o c 'i> with c- e A n _i(Z). If Z is smooth, then the Chern-Fulton classes 
coincide with the Chern classes of the tangent bundle. Since c{T V k) = (1 + H) k+1 
where H G Ak-i{P k ) is the hyperplane class, the degrees of the Segre classes and 
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those of the Chern- Fulton classes are related by 

deg(c:) = E("+ 1 ) deg( Sp ). 

Recall that in the smooth case, the degree of the top Chern class of the tangent 
bundle is equal to the topological Euler characteristic. Thus, in case Z is smooth, 
deg(c' n ) is the topological Euler characteristic of Z and Procedure [T] provides a way 
of computing this topological invariant. 

5. Examples 

In this section we illustrate Procedure [T] with some examples. In these ex- 
amples, for an n-dimensional subscheme Z C P fc , we use the notation cr(Z) = 
(deg(s ),...,deg(s„)), where s(Z,P k ) = J27=a s * and s * e Ai-i(Z). 

Example 5.1. Let / c C[x, y, z] be the ideal / = (x 2 , y 2 , xy) and let p £ P 2 be the 
degree three zero-scheme defined by /. Then s(p, P 2 ) = s £ Ao(p), and A (p) = Z 
via the degree map. Now let /i,/2 £ I be general elements of degree 2 and put 
J = ih,h)- Then J : I = (x,y) and (J : I) : I = (1). Hence J : I°° = (1) and the 
residual R is empty. Therefore 

so = 2 2 - deg(i?) = 4. 

Note that Sq is not equal to the degree of p (which is equal to 3). 

Example 5.2. Consider a plane curve C C P 2 defined by one element g £ C[x, y, z] 
of degree m. Then it is immediate from Procedure [T] that <r(C) = (in, — m 2 ). In 
case g — xy, we get u(C) = (2, —4). 

Now consider the scheme DCP 2 defined by I = (x 2 y,xy 2 ). The support of D 
is the union of the lines L\ = {x = 0} and L 2 = {y = 0}, but D has an embedded 
point at p = LiDL2- The normal cone CdP 2 has three irreducible components, all 
of dimension 2, and the supports of these components are L\, Li andp, respectively. 
The supports are the so-called distinguished varieties of the intersection defined by 
{x 2 y, xy 2 }. A general element j\ £ 1(3) may be written fi — xy(ax + by), for 
general a, b £ C. Then (/i) : I°° — (fi) : I — (ax + by). Hence the residual R is 
the line {ax + by = 0} and 

deg(s ) = 3 - dcg(i?) = 2. 
For a general fa £ 1(3) we have that (fi, fi) = I and it follows that 
deg(ai) = 3 2 -2-3deg(s ) = -3. 

In summary, a(D) = (2,-3). Observe that the Segre classes detect the embedded 
point p. 

Example 5.3. To illustrate Procedure [1] we show in this example how it works on 
a surface ZC P fc . Let Z be defined by an ideal I which is generated by polynomials 
g Q ,...,g r in C[x , ...,x k }- Let to = max^deg^)} and let fi,...,f k £ I(m) be 
general. Let J 2 = . . . , fk-2), h = (h, • • • , fk-i) and J = (/1, . . . , f k ). Let Wi 
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be the scheme defined by Ji and let Ri be the scheme defined by Ji : Then 
W 2 = Z U R 2 dim(i? 2 ) = 2 (or R 2 is empty), 
Wi = Z U Ri dim(i?i) = 1 (or R 1 is empty), 
W = Z U Ro dim(i? ) = (or R is empty). 

The degrees of the Segre classes so, s%, s 2 of Z in P k are computed as follows: 
deg(s ) = m fc - 2 -deg(i? 2 ), 
deg(si) = m fc ~ 1 - deg(i?i) - (k - l)mdeg(s ), 

deg(s 2 ) = m k - deg(R ) - Q^m 2 deg(s ) - fcmdeg(si). 

6. Implementation and benchmarks 

There is an implementation of Procedure [1] in the symbolic setting using the 
software system Macaulay2 [7] and it is available at |http://www.math.su.se/~jost/| 
It uses the improvement of Rcmark [4.1l The user may choose to be given the degrees 
of the Chern- Fulton classes as output. 

As an alternative to the Grobner basis computations carried out in Macaulay2, 
one can use the regenerative cascade algorithm [9] implemented in the software 
package Bertini [3 ]. The regenerative cascade algorithm uses numerical homotopy 
methods to collect data about solution sets of polynomial equations and this data 
includes the degrees of the residuals that are used in Procedure Q] to compute the 
degrees of the Segre classes. 

Table Q] shows run times on some examples, comparing our implementation 
"segreClass" to two other algorithms. One is the April 2009 version of "CSM" 
which implements Aluffi's algorithm to compute Segre classes, see [T|. The other 
is the routine "euler" from Macaulay2 which computes the topological Euler char- 
acteristic of a smooth projective variety. Observe that the input to "euler" is a 
projective variety, not an ideal. Following Table [TJ we provide a few details about 
each example. Additional details on how to generate the equations for each example 
may be found at |http: / /www.math.su.se/^jost /| 

Table 1. Comparison of run times. On an AMD Athlon 64 Pro- 
cessor, 2.2 GHz, and with 1 GB RAM. The computations marked 
with "-" were terminated after 3 hours. 



Input 


segreClass 


CSM 


euler 


Rational normal curve in P B 


0.5s 


180s 


4s 


Rational normal curve in P 10 






512s 


Grassmannian G(l, 5) C P 14 




2s 




Smooth surface in P 8 defined by minors 


89s 






Abelian surface in P 4 


175s 






Segre embedding of P 2 x P 3 in P 11 




8s 





The defining equations of the rational normal curves in Table [T] are given as 
(2 x 2)-minors of a matrix with variables as entries. The Grassmann manifold is 
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embedded with the Pliicker embedding. The surface in P 8 is defined by the (2 x 2)- 
minors of (4 x 3)-matrix of random linear forms. The ideal of the Abelian surface is 
generated in degrees 5 and 6. The Grassmannian and the Segre product were run 
over Q and the other examples were run over the finite field with 32749 elements. 

Remark 6.1. In connection with the comparison made in Table Q] it should be 
noted that the routine "euler" computes the topological Euler characteristic by 
first computing the Hodge numbers of the variety and then taking an alternating 
sum of them. Thus, "euler" computes interesting information that is not attainable 
from the Segre classes in any obvious way. 



7. Conclusions 

This paper presents an elementary algorithm, based on residual intersection, 
to compute the degrees of Segre classes of a subscheme of projective space. The 
algorithm has been implemented in Macaulay2 [7] and can be implemented in a nu- 
merical setting using numerical homotopy methods. It is particularly well suited to 
the regenerative cascade algorithm [5] implemented in the software package Bertini 
[3]. The table of example run times illustrate the complementary nature of the 
symbolic implementation of the algorithm to previous symbolic algorithms, in par- 
ticular to the algorithm of Aluffi PQ in the general case and to the algorithm "euler" 
found in Macaulay2 when run on smooth projective varieties. The possibility of 
a numeric implementation shows promise for extending the range of problems to 
which the algorithm can be applied. 
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